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Frequency Energy Multiplier Approach to Uniform 
Exponential Stability Analysis of Semi-discrete 
Scheme for a Schrödinger Equation 
under Boundary Feedback 


Bao-Zhu Guo and Fu Zheng 


Abstract—In this paper, we investigate the uniform exponential 
stability of a semi-discrete scheme for a Schrödinger equation 
under boundary feedback stabilizing control in the natural state 
space L? (0, 1). This study is significant since a time domain energy 
multiplier that allows proving the exponential stability of this 
continuous Schrödinger system has not yet found, thus leading 
to a major mathematical challenge to semi-discretization of the 
PDE, an open problem for a long time. Although the powerful 
frequency domain energy multiplier approach has been used in 
proving exponential stability for PDEs since 1980s, its use to 
the uniform exponential stability of the semi-discrete scheme for 
PDEs has not been reported yet. The difficulty associated with 
the uniformity is that due to the parameter of the step size, 
it involves a family of operators in different state spaces that 
need to be considered simultaneously. Based on the Huang-Priiss 
frequency domain criterion for uniform exponential stability of a 
family of Cy-semigroups in Hilbert spaces, we solve this problem 
for the first time by proving the uniform boundedness for all the 
resolvents of these operators on the imaginary axis. The proof 
almost exactly follows the procedure for the exponential stability 
of the continuous counterpart, highlighting the advantage of this 
discretization method. 


Index Terms—Schrédinger equation, boundary damping, fre- 
quency domain multiplier, semi-discretization, uniform exponen- 
tial stability. 


I. INTRODUCTION 


Ontrol systems described by partial differential equations 

(PDEs) is infinite-dimensional. Being such, its controller 
such as the observer-based feedback control is also infinite- 
dimensional. As a result, the discretization finds itself in 
almost all implementations of PDE control. Among many 
discretization methods is the finite-difference method which 
becomes popular due to its simplicity in principle and its 
appeal to engineers. One of the most commonly used dis- 
cretization method is the so-called semi-discrete scheme which 
keeps time continuous while discretizing the spatial variable. 
It has been widely studied in literature. The main advantage 
of the semi-discrete scheme is that it results in an ordinary 
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differential equation system, which control researchers are 
most familiar with. However, it has been acknowledged for a 
long time that the uniform exponential stability with respect to 
the spatial discrete step size cannot be guaranteed for classical 
semi-discrete schemes for PDEs, largely due to presence of 
high frequency spurious components. In addition, some other 
typical important control properties such as uniform observ- 
ability and uniform exact controllability cannot be guaranteed 
either. The reason for this loss is that the spurious modes 
are only weakly damped in the process of semi-discretization. 
A detail account can be found in [24]. For wave equations, 
several remedies such as Tichonoff regularization [7], mixed- 
finite elements [2], [17], high frequency filtering [10], and 
non-uniform meshes [4], have been proposed to circumvent 
this difficulty. Among many these remedies, the numerical 
viscosity damping introduced in [19], [20] is the most popular. 
However, this approach brings a viscosity term artificially 
added into the classical discrete scheme. The coefficients of the 
numerical viscosity damping vary from PDE to PDE. Recently, 
a new natural semi-discrete scheme based on order reduction 
finite difference method was introduced in [13] and has been 
applied to different systems [8], [23]. This approach has the 
critical advantages that it guarantees the uniform exponential 
stability. In addition, as a natural semi-discrete scheme, it al- 
lows one to prove the uniform exponential stability in a manner 
parallel to its continuous PDE counterpart. Nevertheless, all 
the previous papers on this scheme involved construction of 
Lyapunov functional which the proof heavenly relies on, both 
for semi-discrete scheme and the continuous counterpart. 

Construction of a suitable Lyapunov functional for a PDE 
relies on a time domain energy multiplier, which is not always 
available and its construction is most often very technical. In 
1980s, a frequency domain energy multiplier approach was 
developed for the exponential stability initially for a single 
PDE ([15]). The approach is based on a frequency domain 
characterization for exponential stability of Co-semigroup in 
Hilbert space. Originally developed independently in [9] and 
[18], the result of was proved later in [16] to be valid for 
uniform exponential stability of a family of Co-semigroups in 
Hilbert spaces as well. Uniform admissibility and observability 
for the finite element space semi-discretizations of abstract 
Schrödinger system and second order infinite dimensional 
vibrating systems have also been developed [5], [6]. 

In this paper, we investigate the uniform exponential sta- 
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bility of an order reduction semi-discrete scheme for a 
Schrödinger equation under boundary control by the frequency 
domain multiplier approach. It is significant because one 
cannot find a suitable time domain Lyapunov functional both 
for the continuous PDE and for its discrete scheme. This 
implies that successful approaches presented in [14], [13], [8], 
[23] cannot be applied here. As a matter of fact, in order to 
apply the Lyapunov method, the paper [14] has to consider 
the Schrodinger system in the high order state space H'(0, 1), 
whereas our state space is the standard space L7(0,1). The 
problem in L7(0,1) has been open for quite a long time. Fairly 
speaking, this paper brings a new way to the proof of the 
uniform exponential stability of the semi-discrete scheme for 
PDEs. It is also worthy pointing out that the proofs for both 
continuous PDE and for the discrete counterpart are again 
analogous, demonstrating the advantage of the order reduction 
semi-discretization approach. 

We proceed as follows. In the next section, Section II, we 
prove the exponential stability of the continuous PDE by the 
frequency domain multiplier method. Although it is the sim- 
plest PDE ever studied in the literature, it helps in constructing 
a frequency domain multiplier for its semi-discrete counter- 
part. In Section III, we design a semi-discretized scheme and 
obtain a family of finite-dimensional systems. In Section IV, 
the uniform exponential stability is developed by the frequency 
domain multiplier approach. We introduce the shadow ele- 
ment to help understand the numerical approximating scheme, 
which plays an important role in the proof of uniform stability. 
Some concluding remarks are included in Section V. 


II. STABILITY OF SCHRÖDINGER SYSTEM VIA FREQUENCY 
DOMAIN MULTIPLIER 


Consider the following Schrödinger equation under bound- 
ary control: 


w;(x,t) = —iwy,(x,t), t > 0, x € (0,1), 
w(0,t) =0, t>0, 


w,(1,t) =u(t), k>0, t>0, (1) 
y(t) =w(1,t), t 2 0, 
w(x,0) = w(x), x € [0,1], 


where u(-) is the control, y(-) is the measured output and wo(-) 
is the initial state. Under the proportional feedback control: 


u(t) = —kiy(t), k>0, (2) 
the closed-loop system of (1) becomes 


w;(x,t) = —iwyx(x,t), t > 0, x € (0,1), 
w(0,t) =0, t>0, 

w,(1,t) = —kiw(1,t), k>0, t>0, 
w(x,0) = w(x), x € [0,1], 


(3) 


We consider system (3) in the natural state space L7(0,1). 
Define the system operator of (3) as follows: 


Af =—if" Vf € D(A), 
D(A) = {f € L° (0,1)|f € H?(0,1), (4) 
f(0) =0, F0) = -ikf(1)}. 
Then, (3) can be written as an evolution equation in L7(0, 1): 
w(-,t) =Aw(.,t), 
oe = K a 


It is seen that 


1 
Re(Af, Pon =Re | -kf Fd =k FDP, © 


which implies that A is dissipative. In addition, the operator 
A is invertible and 


—kx fo xf (x)dx 
1+ki 


-i f &- ofar -i [xp (ax 


which is bounded in Z? (0,1). As a result, A generates a Co- 
semigroup of contractions on L?(0,1) by the Lumer-Phillips 
theorem ([21, Theorem 3.8.4]) and since A~! is compact, the 
spectrum of A consists of isolated eigenvalues only. 
Furthermore, define the system energy for (3) as 


A`' f(x) = 
(7) 


1 fl 
E(t) = 5 | wapat, (8) 
2 Jo 
which is non-increasing as a consequence of (6): 
E(t) =—Kw(1,1)/?. (9) 


We point out that a different version of (3): 
wy, (x,t) = —iwyx (x,t), t >0, x € (0,1), 
w(0,t) =0, t>0, 
wx(1,t) = —kw;(1,t), k>0, t>0, 
w(x,0) = w(x), x € [0,1], 


(10) 


was investigated in [14], for which one can find a time 
domain energy multiplier, and a Lyapunov functional was 
then constructed to both system (10) and its semi-discrete 
counterpart. However, system (3) is a rather unusual system 
for which a time domain energy multiplier has been not 
found yet. A first exponential stability result of system (3) 
was proved by the Riesz basis approach in [12]. Although 
the Riesz basis is powerful and the result obtained is much 
deeper than the result obtained from the multiplier method; for 
instance the spectrum-determined growth condition is usually 
a consequence of the Riesz basis approach yet this is usually 
not the case with the multiplier method. Unfortunately, the 
Riesz basis is extremely difficult, at least at the moment, to be 
applied for the uniformly exponential stability of semi-discrete 
model for (3) developed in this paper. 

In this paper, we use an alternative powerful method called 
the frequency energy multiplier method, which has been 
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developed for continuous PDEs over the last three decades 
({15]). In stability analysis, we can almost give one-to-one 
correspondence from continuous system to its discrete coun- 
terpart using this method. Our approach is so powerful that 
can be applied to other PDEs as well. For notation simplicity, 
hereafter, we omit without confusion the obvious dependency 
in time and spatial domains. The C” denotes the n-dimensional 
complex Euclidean space; the N* stands for the set of the 
positive integer numbers; and R the set of real numbers. 
Since A generates a Co-semigroup of contractions on 
L?(0,1), a well-known result of Hung-Priiss theorem [9], [18] 
states that the Co-semigroup generated by A is exponentially 
stable if and only if it possesses the following two properties: 


1) Every imaginary number belongs to the resolvent set of 
A, that is, iR C p(A). 

2) The inverse operator of i@—A is uniformly bounded for 
all imaginary numbers, that is, 


sup || (i@ =A)" || < æ. 
OER 


(11) 


The property iR C p(A) is stated in the following Lemma 
2.1. 

Lemma 2.1; Let A be defined by (4). Then, iR C p(A). 

Proof. If there exist B € R,B 40 and a nonzero f € D(A) 
such that iB f = Af, then 


iB f (x) =—if" (x), 
f'(1) =—kif(1), f(0) =0. 


Take the inner product with f(-) over [0,1] on both sides of 
the first equation of (12) to obtain 


ipl si? =e eye +i f Pas 


which gives f(1) =0 and hence f’(1) =0. This shows that 
(12) has only zero solution, a contradiction. | 


Theorem 2.1: Let A be defined by (4). Then, (11) holds 
true. As a consequence, the Co-semigroup e^ generated by A 
is exponentially stable in L7(0, 1). 

Proof. We prove by assuming contrary of (11) that there 
exit a sequence @, > ©, fn € D(A), || ful] = 1 that 


(12) 


(13) 


Jim || @@n — A) fall =0 
i.e., 
i@nfa if, — 0 in L’ (0,1). (14) 
Since 
Re((i@n —A) fn, fn) 22(0,1) = Re(—Afns fn) = Kl fn(L)? + 0 
(15) 
by the boundary condition f/(1) = —ikf,(1), it gives 
fl) 30 (16) 


From (14) and ||f,|| = 1, it follows that AG) is bounded in 
L?(0,1). By 


< [lfa ll; 


RO — (1) =| fi Ras 


it follows from (16) and @, — œ that 
1 

Ja) aboni L (0,1). 
On 


Since 


n X 


n a 
Re(Ole tet) = 


1 1 l 
+g ODP- 5 f, Mean, 


and 


we have by (15) and (16) that 
1 1 1 
f h at S OPa, a8 
0 Qn JO 
which shows that when @, > 0, ||fal|? > 0, which is a 


contradiction to || f,|| = 1. On the other hand, since from (14) 
and @, — œ, we have 


f aE zuf 
=f (n+ A) (r FD) a 


= f' [op AOE) a 


HES MOROROOO 
= | | ae OL] ax 


Hno AORTADA] 


dx 


(19) 


I 2 
= dx— 0, 
Zy Lins) Pa 


Substitute f/(1) = —ikf,(1) and f,(0) 
(15)-(16) to obtain 


1 1 f! 2 
[imcrPar+ f Mee 


which shows that when @, < 0, || fal] 3 0, which is also a 
contradiction. | 


= 0 into (19), and use 


2 1 
dx- = | fi (x)2dx 3 0. (20) 
Qn JO 


I? 


II. SEMI-DISCRETE SCHEME OF SCHRÖDINGER EQUATION 


In this section we apply the order reduction method to derive 
a semi-discrete scheme for (3). To this purpose, we introduce 
an intermediate variable v(x,t) = wx(x,t) to reduce the order 
of the spacial derivative of (3). In this way, the Schrödinger 
equation (3) can be rewritten as the following equivalent form: 


wi (x,t) +ivx(x,t) =0, 
v(x, t) = wx(x,t), 
w(0,t) =0, 

v(1,t) = —kiw(1,t), 
w(x,0) = w? (x). 


(21) 
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The semi-discretization process is similar to [14]. For the 
sake of completeness, we sketch briefly the process. For fixed 
N € N*, consider an equidistant partition of interval {0, 1]: 


O=x9 <x] <+ <x = jh <i) < x41 l, 


where h = wt is the mesh size. Denote the sequence {u;}y +! 
by {u;}; and introduce respectively the average operator and 


the first-order finite difference operator as 


Aj j _ #4j+1— Uj 
= up E D 
For the solutions v(x,t) and w(x,t) of (21), let {V;(t)}j and 
{W;(t)}; be grid functions at grids {x;}j, satisfying 

V(t) = v(@j,t), 


W,(t)=w(x;,t), O<f<NH+1. 


The first equation of system (21) holds at (Xj ,t), i.e., 


w(x jp 1st) + ive (eit) =0, 


— + 1 7 a 66999 _ 
where Xl = (j+ z)h. Hereafter the prime “/’ represents 
the derivative with respect to time t. Replace the differential 
operator 0, with difference operator ô, to get 


y(t) = 0h’). (23) 


Similarly, for the second equation of system (21), it has 


Vist (t)— 8&W, 1) = Ol). (24) 


J+3 
By dropping the infinitesimal terms in (23) and (24), and 
replacing W;(t) and V;(t) by w;(t) and v;(t), respectively, we 
arrive at a semi-discretized finite difference scheme of system 
(21) as follows: 


w’. 1 (t) +iðv, 1 (t) =0, 0<j<N, 
ITZ Iry 
vj) =O U O<j<N, 


(25) 


where vj(t) and w;(t) are grid functions at grids x; (0 < j < 
N+1), and wo is the approximation of the initial value w°(x;). 

Remark 3.1: The semi-discretized system (25) is a family 
of differentiation-algebra systems, which is called singular 
systems for which there are huge amount of references related 
to them. See for instance [3], [11], [22] and the references 
therein. 

Now, we eliminate the intermediate variables v;(t) from 
(25). To this purpose, let 


Wi(t) = (wi(t), walt), +++ wwa (E) 
be unknown variable of (25) and 


V(t) = (vot), vi (t); 


the auxiliary variable. We write (25) into vectorial form: 


0 
Dp W; (t) = —iMpV;,(t) — i ; 
khwn4i(t) 
0 (26) 
D; Va(t) = —My Wa(t) + i ; 
i2-'kwy41(t) 
W,,(0) = (wp, we, a why)" 
where the matrices Dp and Mp are given by 
1 
1 1 
1 
Dy, = 2 Eo ’ 
1 
1 (N+1)x(N+1) 
(27) 
—1 1 
-1 1 
M 1 
_ h 5 
-1 1 
-1 (N+1)x(N+1) 


Obviously, both D; and Mp are invertible. The differential 
algebraic system (25) or (26) can be written as an evolution 
equation in CNt!; 


Wi (t) = AWalt), Walt) € Yn = CF, (28) 
W, (0) = (w, w9,- Wea)” E Yn, 
where %, is defined by 
AY p 
af —1 wane 
=D," [iM, (Dj) ' (M7 Ya (0,--- ,0,2-likyny1)") 
- D; ' (0, a ,0,kh-!yw41)"] > 
WYn = (v1, 9250+ yng)! E€ CN, 
(29) 


System (28) is naturally discussed in the state space C'+!, 
To relate C+! in (28) with the step size, we write Y} = C+! 
and define a new inner product for Y}: 


(Yna) =h ( DiYn, Dia) ,VYn, Yn & Yn, 
h 


where (-,-) is the standard inner product of C+!, For Yp = 
(y1,¥2,°:: ,¥wi1)! € Ya, we choose the vector 


Zn = (Z0;Z15°°° zN) ect satisfying 


Dj Zn = —M; Yn + (0, ++: 0,27 Niky)! 
We call Zp the shadow element of Y;,, which can simplify 
significantly the notation in the later proofs. 


The classical semi-discrete scheme is similar with (28) 
where the average operator Dp = Iy+1, 1.€., 


Wi (t) = AW, (t), Walt) € Yr =C, 


Wa (0) = (w, w9, mae Mya)! € Ya, 


(30) 


(31) 
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in which the Â, is defined by Proof. For Y, = (y1,2,--:,yw41) € Yn, let Zp = 
` (z0,Z1,*+* ,ZN) be the shadow element of Yp: 
AnYn = iMn (M; Ya — (0, ,0,2-Nikyw 1)" ) T T 1: T 
(32) D, Zn = —M, Yn + (0,--- ,0,27 ikyn+1) , 
—(0,---,0,kh-"yw41)'. -1;_. -1 T (4) 
AY, = D; [—iM),Zn + (0,--- ,0,kh™ yn +1)’ J. 

At the end of this section we explain the significance of i sa . z 
the discrete scheme (28). We plot two figures in Figures 1 me 20 ae ane A +i h To a8 tas Yn, = 
and 2, respectively. Figure 1 depicts the maximal real parts 01:923: Yi) € Ya such that 4Y = Yn. Then, 
of the eigenvalues of the classical semi-discrete scheme (31) D} Zn + (0, ,0,27!zy41)! = —M}} Yn, (35) 
with step size h, from which we see that the real parts of the P f 
eigenvalues approach zero. Figure 2 depicts the maximal real which is equivalent to 
parts of the eigenvalues of the order reduction semi-discrete zai = 8y, j=0,1, N, (36) 
scheme (28) with the same step size, from which we see that PEZ IVD 
the real parts of the eigenvalues approach a negative number. and 
In both figures, we take k = 1. DY, = —iMpZn — (0,--» ,0,ih!zy41)", (37) 

which is equivalent to 
0 T pe eee FOS Ft 6 A, 
è (38) 


Real part of eigenvalue 


-0.1 b 


Fig. 1. Maximal real parts of eigenvalues of the semi-discrete scheme by 
classical method (31) 


y k k I I e a e E 


x 
+^ 


Real part of eigenvalue 


Fig. 2. Maximal real parts of eigenvalues of the semi-discrete scheme by 
order reduction method (28) 


IV. PROOF OF UNIFORM EXPONENTIAL STABILITY 
This section is devoted to the proof of the uniform expo- 
nential stability of (28). To begin with, we first show that 2%, 


is dissipative for every step size h. 
Lemma 4.1: For the matrix .%, defined by (29), there holds 
Re (AYnsYa)y, =—klywsil?s Y ¥en, (83) 


which implies that æ% is dissipative for every h € (0,1). 


Yih =I ee iid j=0,1,---,N, 


where in all (35) to (38), it was assumed that yo = 0. Take the 
inner product between .%,Y;, and Y, in Y} by taking (36) and 
(38) into account to obtain 

Re (An Yn) y, = Re (Fa Ya) 


Yn 


= h = 
( Din Din) +5 ( DaYn Daa) 


h&n = , 
=5 Ija ty Yj (using (38)) 
j=0 j=0  * 
hi & 2 hi X E (39) 
Sog Moen ty Lopez 
j=0 j=0 
hi & = 
T ) [Bea (0) +2), 18544 
F a 
hi = 2 = 
t5 = [y xZ jl TORRI (using (36)) 
= 
A simple calculation shows that 
hi X T _ 
= | Ting ERETON 
I= 
se ee 
+7 È bis nd + Buy Eja] 
= 
i N 
a y [(zy41 =z) Oja HY) + (Zi +z) Oj -y,)| 
j=0 
- N 
i 2 _ B 7 
t7 Y [Op +99) Em z) + O —y) Zu) +z )] 
j=0 
ið eu 
a D j =a] +5 D Diza —y;zj] 
J= {= 
Vo _ i _ _ 
= 5 [zoo — zw+1¥n4i] + 5 b+ 12Zw+1 — yo%0] 


= —k|yn+1 |7. (using —ikyy+1 =zv+1 and yo = 0) (40) 
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The (39) and (40) leads to (33). | 
Define the energy of (28) as 


N 
Fu) = 5 È | je = 3 MOWO, 


(41) 
which is the discretization of the continuous energy (8). The 
following Lemma 4.2 is the discrete counterpart of (9), which 
is a consequence of (33). 

Lemma 4.2: The E;(t) defined by (41) satisfies 

E,(t) = kwny (O)|. (42) 

The dissipativity of .%, implies that the spectral set o(.%,) 
of æ%, is contained in the closed left half-plane of the complex 
plane C. Actually, we have more stronger result. Precisely, for 
any 0 < h < 1, the spectral set o(.%,) of æ% is contained in 
the open left half-plane of C. This is the following Lemma 
4.3. 

Lemma 4.3: For every h € (0,1), iR C p(%). 

Proof. If there exist B € R and nonzero Y, € Y, such that 
iBYn = GYp, then it follows from (33) that 


0 =Re (ibYn, Yn)y, =Re (A Yn, Yay, =—Alywsil?. (43) 
Replacing Yh by ißYn in (37), we obtain 


By p41 + Oz 541 =O, O<j<QN, 
(44) 
Zj+4 ZEY =0, 


where Zp is the shadow element of Y, defined in (34), yo = 
0 and zy41 := —ikyņy+1. Hence zy+1 = yy+1 = 0 from (43). 
Setting j =N in (44) yields 


Bhyn = 2zy, ZN = -Žyy. 
It follows that yy = zy = 0 whenever Bh? +4 is nonzero. Under 
the condition Bh? +4 40, suppose z j+1 =yj4+1 = 0 and solve 
(44) to arrive at z; = yj =0. This gives Y, = 0 by induction, 
which is a contradiction. On the other hand, whenever Bh? + 
4=0, it follows from (44) that 


1 1 ; 
=a +y) = = (Zj41 z) j=0,1,---,N, 
h 2 
j 1 (45) 
=j — yj) = (Zoi +z), j =0,1,---,N, 
h 2 
which implies that 
h h 5 
Yj = zZj+b YI Si j=0,1,---,N. (46) 


This, combining with yọ = 0 and yy+; =Q, gives yj =0 (j = 
1,2,--- ,N) which is also a contradiction. This completes the 
proof of the lemma. | 


The following lemma comes from [13]. 


Lemma 4.4: Let {uj}i, {vi}; and {w;}; be the sequences of 
complex numbers. Then, 
1 
Z YE (uia — ui) (vii +1) (Witt + wi) 
i=0 


1 N 
Ea (uit1 — ui) (vi+1 — Vi) (Wit — Wi) 
i=0 
Ie (47) 
+3 (ui+1 + ui) (visi — vi) (wi+1 + wi) 
i=0 
1 N 
+2 (uit1 +ui)(vi+i + Vi) (Wit — wi) 


= UN+1VN+1WN+1 — UOVOWO. 


The following uniformly stability criterion which was pre- 
sented in [15] or [1] will be used in the proof of our main 
result Theorem 4.2 later. 

Theorem 4.1: Let h* >0 and let {Si(t)}ne(ons) be a family 
of semigroups of contractions on the Hilbert space Hp, and let 
An be the corresponding infinitesimal generators. The family 
{S,(t)} is uniformly exponentially stable if and only if the 
following two conditions are fulfilled: 

e For every h € (0,h*), iR C p(n); 

° SUP j.c(0,n*) per |I(BI— An) || <o. 

Now, we are in a position to give the main result of this 
paper. 

Theorem 4.2: For the matrices 2%, defined by (29), the 
corresponding family of Co-semigroups Tj,(t) generated by 
Æp is uniformly exponentially stable, that is, there exist two 
constants M > 0 and @ > 0 independent of h € (0,1) such that 


ITRE) < Me“, Vt > 0. (48) 
Proof. The proof is based on Theorem 4.1. Notice that by 
Lemma 4.1, for every h € (0,1), 7),(t) is a Co-semigroup of 
contractions. The fact that .2%, satisfies the first condition of 
Theorem 4.1 has been claimed by Lemma 4.3. In order to show 
that the family .%, satisfies the second condition of Theorem 
4.1, we prove by contradiction. If the second condition of 
Theorem 4.1 is false, then there exist a sequence B, € R, hy € 
(0,1), and Y% € Ya Yp; lly, = 1 such that 
[UR lly, <2) UR, = Balin — Din Yk (49) 
By the Cauchy-Schwartz inequality, it follows from (49) and 
(33) that 


Re (Uh Yh), = Re (Yi Yn) 


Yin Yan 


j (50) 
=Klyn,+il" S n', 
Let 
Zin = (0zi n TA € Zh 
be the shadow element of Yẹ = (Y1,y3 „Yi, +1) (see(30)), 


U} = (ui,u3,-++ uy 44)! with An(Na +1) = 1. Set artificially 
uy = Yo =O and zy.) = —ikyy,., to unify the notation of 
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u” m and êz 
(49) that 


Dj, Uj, _ 
Tyn 
Mhh, 


41 for j =0,1,-:- 
+3 


iPnDr Yh, + iMn,Zh, ok (0, aii 0, ae rey 


T on —lun T 
= Dho Zh, + (0, ,0,2 ZN, +1) , (51) 
or in vector form: 
Wit Yo gi Ôrzo4 1 
n n 
Uil : Jigi ; x2 1 
Z = iBn 2 +i l 2 
n n 6 n 
Z 
Natt Nath Nath 
Zs 1 oxy” 1 
0+3 otz 
n 
zi ôy a 
2 SEN 2. 
z” ô n 
Nath YN, 
(52) 


The proof will be split into three claims and each claim 
corresponds to that in the proof of stability of PDE. Clam 1 
corresponds to @, — œ in the proof of Theorem 2.1. 

Cliam 1: |f,| > C' > 0 for some constant C’ independent 
of nEN*. 

Suppose by contrary that the sequence {B,} contains a 
subsequence which is still denoted by {f,} itself without 


loss of generality converging to zero. Since ||Y} ||y,, = 1 and 

U} lly, <n |, it follows from (52) that 

Nn 2 
h n n at n 
dh Oxi 41 1 Bry 
N, N, 
n n (53) 
n 
< 2h, L w, ” 2B? m Ds 


n —2 
= 2\|U5, ln, + 2Bn Yp IR, <2n", 


which holds for all sufficiently large n. On the other hand, by 
some simple operations, we get 


2 2 
ej zma l = lj j pe — ee ZN 
Nn 2 Nn \ (M > 
=| (az <( YI L leaz 
l=j l=j l=j 
Nn 2 
< (MDI ¥ az 
l= 
Nn 
=m) OnZ ; j=9,1, „Nn, 
jal 
(54) 


in which A(N, +1) = 1 is used in the last step, and for j = 
0,1, Nn 


Il S lj — Anil + emal S q| Ain 3 lz P + lent: 
= 


This inequality, together with ZN, 41 = 
(53), implies that for each j = 0,1,--- 


—ikyy, 4, and (50)- 
Nn ZP = @(n7"). 


Nn. Then, it follows from Therefore, in light of h,(N, +1) = 1 and the second identity 


of (52), 


(55) 
< hn (Nn T 


<Iin(Na+1)Cn7! = O(n !). 
Thus, the deducing process from (53) to (55) tells us that 


Nn 2 
hy Ôz 1 = Oln?), 
jmo 73 
which implies that 
Nn 2 1 
hy L Zal = O(n ) 
j= 


Nn 2 
mY a mY al, 
2a Y jyh = Ziad 
which means, by (55) that 
Nn 
n _ —1 
wd Oxy 1 O(n’) 


Similarly, repeating the procedure from (53) to (55), for Vie 
we obtain 


Nn 


=n 2 
j=0 


which leads to a contradiction. Thus, the sequence {f,,} cannot 
contain a subsequence converging to zero. We can therefore 
assume that |8,| > C’ >0 for some constant C’ independent 
of nE Nt. 

The second claim is the discrete counterpart of (18) but with 
two extra terms 


2 
n a —1/2 
Yith =0(n ), 


2 
LAA 


hW Nn 


n 
4B, 2 


which play important roles in our proofs. 
Claim 2: the following (56) holds true: 


h3 Nn 


Oy, 


Bz 


1 A 
Xi, ln + Bn z in | | Enn Zh lin +2 
lin he 
ZB, aa hn lAn Zi, ree + qin Ann Ye Ile 42 
N, N, 
n y Ma n 
=h, 2 +e P ae 
alin 5 hk? Nn 
taB, OZA [+4 + — 40 Ory 
n d 


= O(n), (56) 
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where (49) and (52) were used. On the other hand, by the 


second identity of (52), we have 


where ||: || cv, +2 denotes the standard norm of C%+? and 
Gn T n \T T 
Zh, = (Z0)215°** 3 ZN +1) = (z) aa) ; 
T T 
Fi, = Oyana) =N) 
1 1 
1 ; 
= 
"2 11 (57) 
LVS wary) 
-1 1 
1 -1 1 
Ah = h , 
-1 1l (N+2)x(N+1) 


Actually, it follows from (49), (52), and DAAA = ] that 
2 


Nn j 
He d a Pr ta - 
N, n n n 
i Ley HAON jad ae 
Mn 
+B SoBe A 


Applying Lemma 4.4 to the two terms of the right hand side 
of (59) and noticing XN, +1 = 1, x =0, x5 
easy to obtain 


ee a 
741747 = Mn, it is 


Nn 
n n n 
Re (i L ÔF ) 


Nn 


= n Ee iaig Ha Dor, Ly; +4 pan 


Br7hn Lo ÂZ ı is uniformly bounded with respect to Nn 
wee teat = 7 Oi +20} HNO} HH) (60) 
j=0 
Nn 
OZ Vues wit +3 L (a HE Oa 2 O +y) 
ô,z” 4 n a u” f j=0 
i 143 Yit BES I+y h Nn 
Bn : Bn : = Inl- n yal -% Le $ 
bz" i u” i 
Natz Natz Natz wid 
Nn 
By (55) and Claim 1, n g 7" 
y (55) 2Re ha D aiT 
Nn Nn Nn 
j =n n n n 
Ba hn De = hn dL H T Hen Dig 
1 
f ; f : Z 61 
is also uniformly bounded with respect to n € N*. Let xj = jhn 4 Lo i +X) +27) Eja —Z) el 
for j =0,1,---,N,-+1 and consider the following estimates: sg 
; +z ih! Api ap a PEt aa — 2) ) 
on 
h yx (Bry 1 + 6,2" \ td N, 
n jb Pny il xe 54 1 n 2 
j=0 J+53 ITZ JTZ Bn = |zn, 411? —hy > sk fre! ay OxZ 
N zn 2 4 j= 
hp Y uw I By (59)-(61), it follows that 
"he jth jth 
j=0 n Nn 3 Nn 2 
n n 
N, a 2 An , A F 4 [+d 
7j 2 j= j= 
< e Vhat, il) (58) i pM 
I n n n 
+e ail hae De, 
MP x 2 Ba jap) Hl Apn Llaza (62) 
<[h y u” Ba “ 
en Oe a A Aer: eA 
I~ I= j+ 
N = Lm, ( wns) 
n 2 F7 t 1 
AA (Biên bi ) i 


= O(n), 


which proves (56) by (50), (58) and zy,41 = —ikyn, +1. 
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The third claim is perfectly the discrete counterpart of (20). 
Claim 3: The following (63) holds true: 


1 5 2 ES 
2 2 2 
Xi, lEn, F p2’ lAn Z, [Enn _ Binley llovn+2 

N, N, 

n 2h, n 
=h z 
nÈ Di yi a Bn o> i+3 

=6(n), 


(63) 
in which Y,, and A; are defined in (57) and ||- ||cw,+2 
denotes the standard norm of C%+ 

Actually, from (52), we have 


IARA 


n 


(64) 


< 
Nie 


On the other hand, it follows from the second identity of (52) 


that z” , = ôy” , and 
its a iF 


Nn 
n 


2hy 
a Z 
B Fie 


n n = 
z" 110Z 1) + 
a H g) 


Nie 


1 7e = n <I = 
= 5B, LO +97) a i) + Oa FY) Ga +7)! 
n j=0 


Nn 


tog BO HNC Z) + 41-9) Gai +27) 


1 Nn z E 
=i E Faza yz) + 
N i 
1 
~ Ba 


Nn 


1 z = 
g L0) 


n -n_n non 
DN +12Ny-+1 FYN, +1ZN, +1 — ¥020 — V020] = 0, 


where zy, +1 = KYN, 41 and yo =0 were used in the last step. 
Hence 

hy Nn = Zi 2 hy Nn 

Br peg p44 +4 jeg) = BF . (65) 
Plugging (65) into (64) and using Oi, lip <n! and Claim 


1, we arrive at (63). 
Finally, if By >0, we have ||Y} ||, = (n!) from (56), 
which is a contradiction to ||¥;"|ly,, = 1. When bn < 0, 


i in = O(n!) by virtue of (63), which is also a con- 


tradiction. We therefore complete the proof of the theorem. 


V. CONCLUDING REMARKS 


In this paper, the uniform approximation of exponential 
stability of a one-dimensional Schrödinger equation is investi- 
gated. We introduce an order reduction space semi-discretized 
finite difference scheme for approximating uniformly the expo- 
nentially stable closed-loop system. Although the scheme has 
been applied to certain PDEs in previous works, they all share 
a common feature that it is possible to find a suitable Lyapunov 
functional for the closed-loop systems, both for the continuous 
system and its discrete counterpart. However, for the system 
considered in this paper, it is a longstanding problem that in 
the natural state space L7(0, 1), even for the continuous system, 
the time domain energy multiplier has not been found. This 
makes the convergence of the semi-discrete scheme of this 
PDE be open for a long time. This paper is the first work 
that applies the frequency domain multiplier approach to the 
uniformly exponential convergence of semi-discretized PDE 
system. The convergence of the discrete scheme to continuous 
system is not included because it is a standard procedure 
and can be followed analogously from [14] and many other 
papers. Considering it is difficult to find a time domain energy 
multiplier for many other PDEs, the approach presented in this 
paper has significant potentials in applying to other PDEs. 
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